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ABSTRACT 

Juniperus deppeana has numerous disjunct populations that include four taxonomic varieties and 
three forms. All four varieties and two forms of Juniperus deppeana from the southwest United States, 
Mexico and Guatemala were analyzed by sequencing nrDNA (ITS), petN-psbM, trnS-trnG, trnD-trnT, 
trnL-trnF. A Bayesian tree gave support for clades of var. deppeana, NM, var. patoniana, and var. 
gamboana. However, several clades with high support contain mixtures of different varieties and forms. 
A minimum spanning network based on 91 mutational events (MEs) showed that the varieties and forms 
are extremely closely related, differing by only 1 to 2 bp (out of 4411 bp). The taxon with the largest 
diffentiation was var. deppeana, Sacramento Mtns., NM that differed by 4 MEs from the Oak Creek 
canyon, AZ individuals. The lack of variation among J. deppeana taxa may be due to the mixing of 
populations during the Wisconsin glacial maximum (70,000 - 13,000 ybp) when life zones descended 
about 800 m. Published on-line: www.phytologia.org Phytologia 95(2): 161-166 (May 1, 2013). 
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Juniperus deppeana Steudel has trunk 
bark that exfoliates in quadrangular plates, thus 
the common name 'alligator bark' juniper. 
Juniperus deppeana is part of the serrate leaf 
margined species of the western hemisphere 
(Adams, 2011) and is widely distributed in the 
southwestern US, Mexico and northern 
Guatemala (Fig. 1). Putative Juniperus d. f. 
sperryi, once known only from the type locality 
in the Davis Mtns., TX, has been found in 
Arizona and New Mexico (Fig. 1). However, 
the bark characters seem to be controlled by only 
a few genes tuhus furrowed bark may have 
arisen independently in western Mexico and the 
southwestern United States (e.g., J. d. var. 
patoniana and J. d. f. sperryi ). Whether the 
furrowed bark trees of Mexico are related to J. d. 
f. sperryi is not well understood. 

The first systematic treatment of the 
serrate leaf-margined junipers was by Martinez 
(1963) who recognized J. deppeana Steudel. var. 
deppeana (checkered bark, (3)4-5(6) seeds/cone, 
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J. d. var. pachyphlaea (Torrey) Mart, (checkered bark, (l)2-4(5) seeds/cone), J. d. var. robusta Mart, 
(checkered bark, (l)2-3(-6) seeds/cone), J. d. var. zacatecensis Mart, (checkered bark, l-4(-7) 
seeds/cone), J. patoniana Mart, (laced bark, (l)2-3(-6) seeds/cone, and J. gamboana Mart, (checkered 
bark, 1(2) seeds/cone) (Fig. 2). 




J. d. var. deppeana, 10539 , El Chico, MX J . cf, var. gamboana. 6863 , Comitan. MX J. d. var. patoniana. TZ 2752. Dur.. MX 


J. d. var. robusta TZ 2832, Dur., MX J. d. f. elongata 10627. Davis Mtns, TX J. d. f. sperryi 11312, Murids Mtn., AZ 



J. d. f. zacatecensis, 6840. Zac., MX 


Variation in bark patterns in J. deppeana 
varieties and forms. 


Figure 2. Variation in bark exfoliation among J. deppeana varieties and forms. 

Zanoni and Adams (1976, 1979) and Adams, Zanoni and Hogge (1984), using morphology and 
essential oils, generally agreed with Martinez's treatment, except J. patoniana was reduced to J. d. var. 
patoniana (Mart.) Zanoni. Additional studies (Adams and Nguyen, 2005; Adams et al. 2007) have 
further clarified geographical variation in J. deppeana. 

Recently, Adams and Schwarzbach (2006) recognized J. gamboana as J. deppeana var. gamboana 
(Mart.) R. P. Adams and J. deppeana var. zacatecensis as J. deppeana f. zacatecensis (Mart.) R. P. Adams. 
Adams and Schwarzbach (2011) found J. deppeana and var. gamboana to be a clade, sister to J. ashei, J. 
saltillensis and J. zanonii (Fig. 3). However, the other J. deppeana varieties and forms were not included 
in their study. 

The focus of the present study was to examine relationships among all the recognized (Adams 
2011) varieties and formas of J. deppeana (except the very minor variant f. elongata , Adams 2011) using 
data obtained from sequencing of nrDNA (ITS), petN-psbM, trnS-trnG, trnD-trnT and trnL-trnF. 
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MATERIALS AND METHODS 

Specimens used in this study: J. 
deppeana var. deppeana, Adams 10539- 
10541, El Chico National Park, 

Hidalgo, MX; Adams 7632-7634, 

Sacramento Mtns., e of Alamogordo, 

NM, USA; Adams 10640-10642, Oak 
Creek Canyon-Flagstaff, AZ; J. 
deppeana var. gamboana, Adams 6863- 
6867, Comitan, Chiapas, MX; J. 
deppeana var. patoniana, Adams 6836- 
6839, km 152, w. of Durango (city), 

Durango, MX (P); J. deppeana var. 
robusta, Adams 10255-10256, w of La 
Ciudad, Durango, MX; J. deppeana f. 
sperryi, Adams 10626, Bridge Spring, 

Davis Mtns., TX, USA; Adams 11312, 

Munds Mtn., AZ; /. deppeana f. 
zacatecensis, Adams 6840-6842, 18 km 
w. Sombrette, Zacatecas, MX; /. 
virginiana, Adams 10231-10232, 

Knoxville, TN, USA. Voucher 
specimens are deposited at BAYLU 
herbarium, Baylor University. 


One gram (fresh weight) of the foliage was placed in 20 g of activated silica gel and transported 
to the lab, thence stored at -20° C until the DNA was extracted. DNA was extracted from juniper leaves 
by use of a Qiagen mini-plant kit (Qiagen, Valencia, CA) as per manufacturer's instructions. 

Amplifications were performed in 30 pi reactions using 6 ng of genomic DNA, 1.5 units Epi¬ 
centre Fail-Safe Taq polymerase, 15 pi 2x buffer E (petN, trnD-T, trnL-F, trnS-G) or K (nrDNA) (final 
concentration: 50 mM KC1, 50 mM Tris-HCl (pH 8.3), 200 pM each dNTP, plus Epi-Centre proprietary 
enhancers with 1.5 - 3.5 mM MgCl 2 according to the buffer used) 1.8 pM each primer. See Adams, 
Bartel and Price (2009) for the ITS and petN-psbM primers utilized. The primers for trnD-trnT, trnL-trnF 
and trnS-trnG regions have been previously reported (Adams and Kauffmann, 2010). 

The PCR reaction was subjected to purification by agarose gel electrophoresis. In each case, the 
band was excised and purified using a Qiagen QIAquick gel extraction kit (Qiagen, Valencia, CA). The 
gel purified DNA band with the appropriate sequencing primer was sent to McLab Inc. (San Francisco) 
for sequencing. Sequences for both strands were edited and a consensus sequence was produced using 
Chromas, version 2.31 (Technelysium Pty Ltd.) or Sequencher v. 5 (genecodes.com). Sequence datasets 
were analyzed using Geneious v. 5.4 (Drummond et al. 2011), the MAFFT alignment program and the 
PAUP* program, version 4.0b 10 (Swofford 2003) for neighbor joining, parsimony, and maximum 
likelihood tree searches. Further analyses utilized the Bayesian analysis software Mr. Bayes v.3.1 


Serrate leaf margined 
Juniperus of North 
America 

Bayesian Tree 
combined 4351 bp 
Nuclear: nrDNA 
Chloroplast: 
petN-psbM, 
trnD-trnT, 
trnL-trnF, 
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Figure 3. Bayesian tree of the serrate Juniperus of North America. 
Numbers at the branch points are posterior probabilities (as percent). 
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(Ronquist and Huelsenbeck 2003). For phylogenetic analyses, appropriate nucleotide substitution models 
were selected using Model test v3.7 (Posada and Crandall 1998) and Akaike's information criterion. 
Minimum spanning networks were constructed from mutational events (ME) data using PCODNA 
software (Adams et al., 2009; Adams, 1975; Veldman, 1967). 

RESULTS AND DISCUSSION 

Sequencing the five gene regions 
(nrDNA (ITS), petN-psbM, trnS-trnG, 
trnD-trnT, trnL-trnF) resulted in 4411 bp 
of data. A Bayesian tree based on these 
data (Fig. 4) supports clades of var. 
deppeana, NM, var. patoniana, and var. 
gamboana. However, several clades with 
high support contain mixtures of 
different varieties and forms. Note that 
accessions of var. robusta are in different 
clades, as are the samples of f. sperryi 
from AZ and TX. It should be noted that 
hybridization between varieties and 
forms should be expected as these 
dioecious taxa are out-crossing plants in 
populations where several infraspecific 
taxa are often present. Figure 4. Bayesian tree of J. deppeana varieties. Numbers at 

the 

branch points are posterior probabilities (as percent). 

Although the Bayesian tree indicates high 
support for some clades (Fig. 4), the magnitude 
of differentiation among accessions is not 
apparent. A minimum spanning network based 
on all mutational events (MEs) assimilates both 
the nucleotide substitution and indel 
information. This analysis revealed 97 MEs, 
with 6 MEs found only once, and 91 MEs found 
multiple times. Of the 91 MEs, 70 differentiated 
J. virginiana from J. deppeana (Fig. 5). Thus, 
the entire differences among these 4 varieties 
and 2 forms amount to only 21 MEs. In general, 
only 1 or 2 MEs separate individuals (Fig. 5). 

The exception is the differentiation of plants 
from NM (Fig. 4) that are separated from the AZ 
plants by 4 MEs. 

Adams and Schwarzbach (2012) found 
that traditionally recognized taxonomic species 
differed by 8 to 12 (or more) MEs, whereas 

varieties appeared to differ by less than 8 MEs. Figure 4. Minimum spanning network. Numbers 
The differences between var .deppeana from AZ on the lines are the number of MEs. 
and NM, although interesting, do not appear to be 
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sufficient, by themselves, to warrant the recognition of a new variety from NM. Additional research on 
geographic variation in leaf terpenes and morphology is in progress so as to examine differentiation in the 
southwestern US. 

CONCLUSION 


The mixtures of various taxa 
within clades may be due to ancient 
climate and past distributions of J. 
deppeana. Wells (1966), using data 
from rat middens from the Big Bend 
of Trans-Pecos, Texas, concluded 
that during the Wisconsin (70,000 - 
13,000 ybp) life zones descended 
about 800 m leading to the 
formation of a piny on-juniper 
woodland in the present Chihuahuan 
desert between the Big Bend of 
Trans-Pecos, Texas and the city of 
Del Rio. Assuming that the effects 
of glaciation were mediated 
southward into Mexico so that life 
zones descended only a few hundred 
meters in Hidalgo, it appears that 
most of the now disjunct 
populations of J. deppeana may 
have once been connected in a 
nearly continuous population of 
distribution around the Chihuahuan 
desert (Fig. 5). It is likely that desert 
peaks within the area concerned also 
supported stands of J. deppeana. 

Wisconsin populations would have 
become spatially separated as dryer, Figure 5. Possible range of J. deppeana during the Wisconsin 

warmer climate developed during the glacial maximum (based on Wells, 1966). The present day 

Holocene (past 13,000 y). Of course, disjunct populations were likely continuous in the foothills 
the Wisconsin was only the most recent around the Chihuahuan desert during the Wisconsin, 
of several pluvial events during the 

Pleistocene, spanning 1.8 my (Flint, 1971). It is likely that during any one (or several) of these pluvial 
events, Juniperus deppeana occupied lower elevation and more southward habitats, leading to more 
contiguous populations in Mexico and the southwestern United States. If divergent populations (or 
varieties) became sympatric during the Wisconsin, this would have facilitated infra-specific crossing. 
This may account for the large genetic variation within some populations. In addition, the millennia of 
continuous populations could explain the lack of differentiation between the recently (Holocene) 
geographically isolated populations. 


I I J. d var deppeana 
] J. d. var gamboana 
P J, d- var, patoniana 
J. d. var. robusta 
SJ.d.f sperryi 
E J. d. f. elongata 
□ J d. f zacatecensis 

J. deppeana range 
/ (ca.) at Wisconsin 

/ glacial maximum 
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